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ABSTRACT 

Representing the semi-hard partonic interactions by the ex- 
change of Lipatov's perturbative Pomeron, we express the semi- 
hard nuclear cross section as a self shadowing cross section. 
With the help of a generating functional technique, we obtain 
average numbers of multiple semi-hard partonic collisions with- 
out any need of using explicit expressions for the multi-parton 
distributions. The average number of semi-hard interactions of 
a given projectile parton against a target nucleus is estimated 
quantitatively and it is shown to grow very rapidly above one 
with increasing the cm. energy. 







1. Introduction 

One of the main problems in the physics of semi-hard interactions is to take 
into account unitarity corrections. The semi-hard regime, in fact, although pertur- 
bative in the elementary interaction, gives rise to cross sections that may violate 
the unitarity limit. Directly related features are the large size of the QCD-parton 
model cross section to produce large p t partons, which is obtained when the cut-off 
p™ lTl , defining the perturbative regime, is lowered down to the mini-jet region[l], 
and the large value of the mini-jet cross section, which has been observed by 
UA1[2]. The standard way to unitarize the semi-hard cross section is by eikonal- 
izing the elementary semi-hard partonic inter action [3]: The semi-hard component 
of the hadronic interaction is included in the eikonal phase as an additive term, 
which represents the inclusive cross section to produce large p t partons, integrated 
down to the lower limit p™ %n and depending explicitly on the impact parameter of 
the hadronic collision. The cross section, including all inelastic events with at least 
one semi-hard partonic interaction, may be easily derived and expressed as a series 
of multiple semi-hard partonic interactions. The resulting distribution in multiple 
semi-hard partonic interactions is a Poissonian, with average number depending 
on the impact parameter and representing different pairs of partons interacting 
independently at different points of the transverse plane. The underlying physical 
assumption is a Poisson distribution for the input multi-parton distributions to the 
multi-parton interaction and, in addition, the absence of any semi-hard rescatter- 
ing for each parton which has interacted already once with semi-hard momentum 
transfer exchange [4]. 

To include semi-hard interactions in high energy nuclear collisions, the eikon- 
alized expression of the nucleon-nucleon interaction can be used as a input to 
construct the nuclear interaction[5] according with Glauber and Gribov[6]. With 
respect to the input parton distributions to the semi-hard component of the nu- 
clear interaction, both possibilities of input nuclear parton distributions, with 
and without shadowing corrections, have been considered [5]. Alternatively [7], the 
whole semi-hard nuclear interaction has been expressed by a Poisson distribution 
of multiple semi-hard parton collisions, as in the hadron-hadron case. The only 
difference, with respect to the case of semi-hard hadron-hadron interaction, is in 
the function of the impact parameter, which describes the overlap of the matter 
distribution in the transverse plane, and in the nuclear parton distributions, which, 
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for nucleus-nucleus, include shadowing corrections [8]. In both of these approaches 
a parton is allowed to interact, with momentum transfer larger than only 
once. Semi-hard parton rescatterings, namely multiple collisions of a given par- 
ton, each with momentum transfer larger than p™ in ? are typically estimated to be 
a negligible effect up to energies of order of lOTeV in the nucleon-nucleon cm. 
system [8]. In this respect the very different conclusion, that semi- hard parton 
rescatterings are, on the contrary, one of the major features in nuclear collisions 
already at energies of the order of ITeV in the nucleon-nucleon cm. system, has 
been drawn [9]. This different conclusions is the origin of the rather peculiar physi- 
cal picture of high energy nuclear interactions where most of the transverse energy 
produced is the result of multiple semi-hard parton collisions [10]. 

The purpose of the present note is to gain a better insight into this last point; 
namely the amount of semi-hard parton rescatterings in very high energy nuclear 
collisions, which appears as a critical element in the physics of nuclear interactions 
at very high energy. We will in fact review the arguments of ref.[9, 10] making 
clear which are the underlying assumptions and, by means of an approach which 
is much more general with respect to the approach used there, we will produce 
estimates of the average number of semi-hard rescatterings which are of much 
wider validity. In fact our result will not depend on any specific form for the input 
multi-parton distributions. 

The paper is organized in two parts: In the first we comment on the general 
framework to discus the semi-hard cross section and we derive the quantities of 
interest for the actual discussion. In the second part we present a few numerical 
results and our conclusions. 

2. General framework and average number of rescatterings 

The inelastic channels are obtained from the nuclear cross section, expressed 
according with Glauber and Gribov, by the use of the cutting rules by Abramovskii, 
Gribov and Kanchelifll]. The rules connect different cut diagrams, which are 
shown to be proportional to one another, in such a way that the evaluation of the 
contribution of the various inelastic channels to the inelastic cross section simpli- 
fies enormously. Actually the result, which is obtained after summing all different 
cut diagrams, is that the cross section acquires the simple probabilistic meaning 



2 



of sum of multiple inelastic interaction probabilities [12]. Remarkably the proba- 
bilistic picture is the quantum mechanical result of the sum of a very complicated 
set of cut amplitudes. The systematic use of the cutting rules allows also to recog- 
nize the validity of general properties of the interaction, such as the possibility to 
single out the self shadowing cross sections[13]. For these quantities the unitarity 
corrections have the peculiar propriety that all contributions from all processes, 
different from the one which is considered, cancel. A stronger statement is that 
also the functional form of the unitarity corrections for a self shadowing cross 
section is given, once the forward nuclear amplitude is expressed in terms of "ele- 
mentary" nucleon-nucleon amplitudes [14]. The question of interest for the present 
discussion is whether semi-hard interactions are quantities of this kind. Two are 
the requirements which are to be satisfied. The first is that the soft corrections to 
the semi-hard process have to cancel, when all contributions, real and virtual, are 
taken into account. The proof of factorization at the level of power corrections [15] 
shows that this is the case when inclusive processes are considered. The second 
requirement, which allows the probabilistic picture of the interaction, is that the 
cutting rules, relating different cut diagrams, are to be valid also for semi-hard 
interactions. While this proof is missing in the general case, it has nevertheless 
been obtained for one of the components of the semi-hard cross section which is 
leading in the high energy limit [16]. One may therefore argue that it is a good 
working hypothesis to assume that, in high energy nuclear collisions, semi-hard 
interactions self shadow, in such a way that the semi-hard nuclear interaction 
can be expressed in a probabilistic way as a function of "elementary" semi-hard 
parton-parton interaction probabilities. The physical picture of the semi-hard nu- 
clear interaction discussed in ref.[9] and in ref.[10] and the estimates of the rates of 
semi-hard parton rescattering, which we are going to present here, are done under 
this basic hypothesis. 

The self-shadowing hypothesis and the probabilistic picture which follows are 
not sufficient to write the semi-hard nuclear cross section. In fact, differently with 
respect to soft interactions, where the number of " elementary" interacting objects, 
which in that case are the nucleons, is fixed, when considering semi-hard interac- 
tions and partons as elementary interacting objects, the number is, on the contrary, 
varying. Actually the configurations with many interacting partons involve the 
multi-parton distributions [17] that are non perturbative quantities, independent 
of the single-parton distributions of large p t physics. The complete description of 
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the nuclear semi- hard interaction needs therefore, as a input, the infinite set of 
non perturbative quantities, represented by the multi-parton distributions. Nev- 
ertheless, for the more limited program to evaluate a few average quantities, much 
less information is needed as a input. In this case it is very convenient to approach 
the problem by means of a generating functional formalism which allows to avoid 
explicit expressions for the multiparton distributions[18]. 

To construct the generating functional for the multi-parton distributions one 
may start with the exclusive multi-parton distributions, namely, at a given scale 
provided by the cut off the probabilities for each of the configurations with a 
given number of partons. The {u\ . . . u n ) are then the exclusive n-body parton 
distributions, where Ui = {bi.Xi) represents the transverse partonic coordinate {pi) 
and longitudinal fractional momentum (a^). The generating functional Z[J] is 
defined as: 

Z[J] = Y J \ I JM...J(u n )W^(u 1 ...u n )du 1 ...du n , (1) 

n U - J 

with the normalization condition 



Z[l] = 1. (2) 

While the exclusive distributions are obtained from the expansion of the generating 
functional around J = 0, the inclusive distributions are obtained by expanding 
Z[J] around J = 1. Explicitly the one-body and two-body inclusive distributions 
and are given by: 



dW(u) =W il \u) + J W^(u,u')du' + ^ J W (3 \u,u',u")du'du" + ... 



5Z 



SJ(u) 



j=i 



D^{ Ul ,u 2 ) =W^(u u u 2 ) + J W^\ Ul ,u 2 ,u')du' 

+ Wi4) (ui,U2,u',u")du'du" ... 



d 2 Z 



5J(ui)5J(u2) 



,/=i 



(3) 
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where the two body correlation (ui, ^2), which measures the deviation from 
the Poisson distribution[18], has been introduced. 

The semi-hard nucleus-nucleus cross section oh is constructed by multiplying 
the exclusive parton distributions of the two nuclei, A and £?, by the probabil- 
ity of interaction, which, after the assumption that semi-hard interactions self- 
shadow, can be constructed from the elementary parton-parton interaction proba- 
bility a(ui, u'j), which represents the probability for the parton i of the A- nucleus 
to have an hard interaction with the parton j of the B-nucleus: 



<j H = J d(3a H ((3) 



J n\ 5J[u\) dJ{u n ) 



x 



^ ml SJ'(u[ -/?)••• 5J'{u' m - P) Zb[J ' ] (4) 

n m 

x i 1 - n n i 1 - u 'j)} } n ^ \ J=JI=0 
i=ij=i 

Here (3 is the impact parameter between the two interacting nuclei and the semi- 
hard cross section is constructed by summing over all possible partonic configura- 
tions of the two interacting nuclei (the sums over n and m) and, for each configura- 
tion with n partons from A and m partons from B, summing over all possible multi- 
ple partonic interactions. This last sum is constructed by asking for the probability 
of no interaction between the two configurations (actually Y\7=i YYJLi^ ~ ^m'] )• 
The difference from one of the probability of no interaction gives the sum over all 
semi-hard interactions. To obtain the average number of semi-hard partonic col- 
lisions one expands the interaction probability as a sum of multiple interactions. 
Since the interaction probability is multiplied by a symmetric expression, one can 
make the replacement : 



n m Q 

i-nnt 1 -^]^^ 1 - n*"i 

t=lj=l N=l , v 

Q (5) 



= S J2 ( n )*i---*nO--&n+i)---0--&q) 

N=l ^ ' 
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where the index N counts the possible interactions, in such a way that Q = nm, 
S is a symmetrizing operator [9], and, in the second line in Eq.(5), the interaction 
probability has been expressed as a sum of multiple interactions. The average 
number of partonic collisions at a given value of Q, (N)q, is given by: 

(N) Q = S^N(Qy i ...* N (l- a N+1 ) . . . (1 - a Q ) (g) 
= SQai 

When in the expression for cr#(/3), as given in Eq.(4), the interaction probability 
is replaced with (N)q, the overall average number of parton interactions (N((3)) 
is obtained. Since SQ&i = mn&i t i, as a result of the the sums on m and n both 
arguments of Za and of Zb are shifted by one unit. One can then write: 



6 J(u{) 



z A [J+i] 



j=o 5J'(u' l -f3) 



Zb[J' + 1] &(m, u'-^dmdu'-L 



J'=0 



(7) 

The overall average number of partonic collisions is therefore expressed, on quite 
general grounds, by the single scattering term, where the parton structure of each 
interacting nucleus enters only with the inclusive one-body parton distribution: 



(N(0)) = J D { ) L \m)a{u u u' 1 )D^\u' 1 -f3)dmdu' 1 = D%\u 1 )®a 1 , 1 ®D%\u' 1 -0) 

(8) 

An analogously general result can be obtained for the average number of 
collisions of each parton. The sum over m in Eq.(4) can be performed explicitly 
in such a way that cth{0) is expressed as: 

'* w = /?5i3S)-u5) w 

n n (9) 

In Eq.(9) every configuration with a given number of partons from the A nucleus 
is required to interact at least once with nucleus B. The probability for one of the 
partons form A, which we label with i, not to interact with the whole nucleus B is 
Zb [l — vim, •)] . One minus this probability gives the probability for that parton 
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to interact at least once. Let us expand this interaction probability as a sum on 
successive interactions: 



m / \ 

X ^ ( i ) • • • ^(1 ~ °i,k+l) •••(!- £*,m) U du' 
k=l ^ ' 



the average number of interactions of the parton under consideration is obtained 
by multiplying by k each term in the sum on k. Both sums on k and on m can 
be done explicitly in a way similar to the case already considered. The resulting 
average number is expressed as: 

(k( Ui )) = J D%\u[-p)a(u i ,u[)du[ = a i , 1 ®D^\u' 1 -p) (11) 

where the structure of the target nucleus enters only with 

One can notice that, while the average number of interactions of the parton 
with the target does not feel more than the inclusive one-parton distribution of 
the target, if one looks instead to the average number of rescatterings of the same 
parton, one is sensible to the multi-parton correlations. Let us in fact evaluate the 
average number of rescatterings (r), namely the average number of interactions 
when the number of interactions is at least two. For a fixed configuration with m 
target partons one needs to compute: 

{r{Uj))m =y^M , Jfr*,l • • -^,fc(l ~ Vi,k+l) •••(!- <Tj, m ) 

k=2 ^ ' (12) 

=ma ij i - m<Tj 5 i(l - a ij2 ) ... (1 - a^ m ) 

The average number of rescatterings (r for an incoming parton with a given 
kinematical configuration represented by Ui, is obtained by multiplying the ex- 
pression in Eq.(12) by the sum on m, giving all the different multi-parton config- 
urations, and by performing the sum. The actual result is: 

(r(u i ))=a i ^[Z B [J'+l]-Z B [J' + l-a i ,]} (13) 

By expanding Zb[J' + 1 — a] for & small and keeping only the first term different 
from zero one obtains: 
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(r(«0) = ® a,, 2 ® [Dg>(«i - ^)d2 } (u^ - 0) + ^CfVi - A «2 " /*)] (14) 

(21 

in such a way that combining (r(ui)) with (k(ui)) the correlation is obtained. 

A similar observation is valid for (N(f3)): By starting the sum in Eq.(6) from 
N = 2, one identifies the average number (y{(3)), which counts those semi-hard 
collisions with at least two semi-hard parton interactions. At the lowest order in 
a, (v{0)) is expressed as: 



(15) 

In Eq.(15) the dominant term is obtained by neglecting the correlations in the 
expression D {2 \u u u 2 ) = D^( Ul )D^\u 2 ) + l/2C^ 2 \u u u 2 ) and keeping only 
the term with the larger number of D^'s: 

(v(/3)) « D { 2 ] <g> cr hl <g> D^D™ <E> a 2)2 ® (16) 

which represents the independent semi-hard scattering of two uncorrelated parton 
pairs localized at two different points in the transverse plane, as it can be seen 
by looking at the dependence on the transverse coordinates in Eq.(16). Actually 
this is the expression that has been used first to make predictions on the rate 
for double parton collisions, in high energy hadronic interactions[19], and later to 
analyze the actual experimental signal[20]. In Eq.(15) the most important term 
containing correlation is: 

D { 1 ] <g> D { 1 ] <g> 5- M ® °2,2 ® + A <-> B 

One may notice that (r(ui)) and give a different information about the 

correlation. In the first case the transverse coordinate of the two correlated partons 
differ only by an amount whose scale is given by the range of a, as it is seen 
by looking at the convolutions in Eq.(14). In the second case the correlation is 
averaged over the whole transverse plane. 
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3. Numerical estimates 

The analysis of the previous paragraph shows that, although multiple parton 
collisions involve a whole infinite set of unknown non-perturbative inputs, namely 
the multi-parton distributions, to construct average quantities, only a limited 
amount of input information is needed. Actually the average number of partonic 
interactions in an inelastic event with given impact parameter j3, (N(f3)), and the 
average number of semi-hard interactions of a projectile parton, with given trans- 
verse coordinate and fractional momentum ttj, against a target nucleus, (k(ui)), 
are both constructed from the one-body inclusive parton distribution 

In order to perform any estimate we need, as a input, When the "ele- 

mentary" interaction probability is small, semi-hard rescatterings can be neglected. 
In this case the integral on f3 of (N(f3)} is equal to the semi-hard cross section 
multiplied by the multiplicity of partons which have suffered a semi-hard interac- 
tion^], and, therefore, it is equal to the inclusive cross section. As a consequence, 
in this case one is allowed to identify with the nuclear structure function 
which is used to describe large pt physics. Notice that this conclusion holds at any 
order in the parton correlations of the multi-parton inclusive distributions. It is 
therefore valid for the particular case of the Poissonian distribution in multiplicity 
of the semi-hard parton collisions, which corresponds to the limit of neglecting 
all correlations[18]. The identification of D^ 1 ' is much less clear when semi-hard 
parton rescatterings are taken into account. If looking to deep inelastic scattering 
on a nucleus, in the kinematical regime of shadowing corrections, the fluctuation 
of the virtual photon in a qq pair lasts long enough to cross the whole nucleus [21]. 
In this case, the interaction of the photon is hadronic, or, more precisely, it can 
be described as the sum of qq states, with frozen transverse distance interacting 
according to the Glauber expression for hadron-nucleus amplitude in the forward 
direction [21]. On the other hand, in the Glauber expression, shadowing is the 
result of all multiple rescatterings of the incoming state, both soft and semi-hard. 
When the amount of semi-hard rescatterings is sizeable, the problem which rises 
is that the quantity entering in the total cross section is the probability to inter- 
act at least once, while the quantity of interest for the present purpose is rather 
the average number of interactions. Nuclear structure functions, with shadowing 
corrections, are therefore an underestimate for On the other hand nuclear 
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structure functions without shadowing corrections are an overestimate, since also 
soft interactions alone produce shadowing at low x. 

To have a quantitative feeling of the importance of semi-hard rescatterings, 
we compute , as given in Eq.(ll), using as a input for both nuclear 

structure functions with and without shadowing corrections. The two choices rep- 
resent a lower and an upper bound for (k(ui)). To perform the actual calculation 
we use the set E and set B structure functions by HMRS[22] and, as scale fac- 
tor, we consider two possible choices, p 7 t ntn /2 and p™ %n . Shadowing corrections 
have been introduced by parametrizing the ratio of the nuclear-nucleon structure 
functions as in ref.8. 

As a further input, we need to specify the elementary interaction. The regime 
of interest is the one of mini-jet production. In the typical configuration two mini- 
jets, namely jets with transverse momentum of a few Gev, have a large separation 
in rapidity. As one expects from asymptotic estimates [23], the rapidity interval 
is filled with radiated gluons, in such a way that higher order corrections in 
play an important role. Asymptotically one obtains that, in the inclusive cross 
section, the first correction to the lowest order result is of order (otsy) 2 [2A], with 
y the rapidity interval between the two interacting partons. To keep into account 
the main features of higher order corrections, we parametrize the "elementary" 
partonic interaction by the exchange of the Lipatov's perturbative Pomeron[23]. 
More explicitly we use the expression of the inclusive cross section to produce 
mini-jets derived by Mueller and Navelet[24]. 

The inclusive cross section to produce mini-jets with transverse momentum 
larger than bGeV, in pp collisions with cm. energy range 200GeV < y/s < 
900GeV, is compared with the cross section measured by UA1[2] in fig.l. Using 
as a input HMRS(B) structure functions and pY lin /2 as scale factor (continuous 
curve in the figure) the gross figures of the data are reproduced. The other choices, 
which we have considered, underestimate the cross section by roughly a factor two. 
The average number of interactions, that a given parton undergoes with the target 
nucleus, (k(b, x)), as expressed by Eq.(ll), is shown in fig. 2 as a function of the 
fractional momentum of the projectile parton x and at zero impact parameter 
b. The two continuous curves refer to the two different choices of the cut off, 
p min = 5GeV and p rmn = QQ e y ^ together with HMRS(B) structure functions 
and p 7 t nin /2 as scale factor, which, according with the UAl data, is the favorite 
choice for the input parameters. The dependence on the scale factor is shown by 
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the dashed line in the same figure, where the input is p™ 171 = 5GeU, HMRS(B) 
structure functions and r p™ in as scale factor. With the same input, the effect of 
using structure functions, with shadowing corrections included, is shown by the 
dotted line. The dependence on the cm. energy y/s is shown in fig. 3, where the 
value of the cut-off has been fixed to p™ 171 = 5GeV. Solid, dashed, dotted lines 
refer to the three values for the fractional momentum of the projectile parton 
x = 1, x = 0.5 and x = 0.25 respectively. In fig. 4 we show the dependence on 
the impact parameter at fixed -/s and p™ m . In all cases the atomic mass number 
is A = 208 and the nucleus is represented as an uniform sphere of radius R with 
sharp boundaries. 

Our conclusion is that the average number of semi-hard collisions, that an 
incoming parton undergoes with a target nucleus, is sizeable larger than one, also 
keeping the lower cut-off p™ TO as high as 5GeV. We mainly limit our analysis to 
the case of interactions between nuclei with sub-energies of the order of ITeV in 
the nucleon-nucleon cm. system. The input to the nuclear case is therefore con- 
strained by the comparison with the available experimental result for the inclusive 
cross section to produce mini-jets, as measured by UA1[2] in pp collisions at similar 
energies. A reason for the big effect, in the number of semi-hard parton rescat- 
terings, is in the large size of the mini-jet cross section which has been observed 
by UAl: almost 20mb at 900GeV cm. energy in pp collisions, with the cut-off 
p-min _ §Q e y Nevertheless an important role is played by the actual represen- 
tation of the "elementary" partonic interaction. When the "elementary" partonic 
interaction is represented by the exchange of Lipatov's perturbative Pomeron, the 
large value of the cross section observed by UAl is the result of a relatively large 
"elementary" cross section rather than of a large flux of incoming partons. In a 
nucleus the probability of semi-hard parton rescattering is therefore enhanced. 
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Figure captions 

1 Inclusive cross section to produce mini-jets with transverse momentum larger 
than 5GeV, in pp collisions as a function of the cm. energy. Continuous line: 
HMRS(B) structure functions, scale factor p 7 t nin /2. Dashed line: HMRS(B), 
scale factor p™ 171 . Dotted line: HMRS(E), scale factor p™ in /2. Experimental 
data from UAl (ref. [2]). 

2 Average number of interactions, (k(b = 0, x)), Eq.(ll) in the text, as a 
function of x in nuclear collisions with sub-energy, in the nucleon-nucleon 
cm. system, y/s = ITeV. Continuous lines: cut off p™ in = 5GeV and 
p min = QGeV, HMRS(B) structure functions, scale factor p™ in /2. Dashed 
line: cut-off p™ in = 5GeV, HMRS(B) structure functions, scale factor p™ in . 
Dotted line: HMRS(B) structure functions, with shadowing corrections in- 
cluded, cut-off pY 1171 = 5GeV : scale factor p rnin /2. Atomic mass number of 
the target A = 208. 

3 Dependence of (k(b = 0,x)) on the nucleon-nucleon cm. energy y/s. Cut-off 
pmin _ §Q e y^ HMRS(B) structure functions, scale factor p rritn /2. Solid line 
x = 1, dashed line x = 0.5, dotted line x = 0.25. Atomic mass A = 208. 

4 Dependence of (k(b,x = 1)) on the impact parameter. Nucleon-nucleon cm. 
energy = ITeV, cut-off pf in = 5GeV, HMRS(B) structure functions, 
scale factor p 7n ' in /2, atomic mass A = 208. 
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